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For intermediate Coulomb energy to Fermi energy ratios rs, spinless fermions in a random 
potential form a new quantum phase which is different from the Fermi glass and the Wigner 
crystal. From a numerical study of small clusters, we show that this phase is characterized 
' by an ordered flow of enhanced persistent currents and occurs for rs values where a metallic 

^\ phase has been observed in two dimensions. 
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1 Introduction 

^ ; 

^ ' A crucial parameter for a system of charged particles is the Coulomb energy to Fermi energy ratio 

• r^. In a disordered two-dimensional system, the ground state is obvious in the limiting cases. 

For large r^, the charges form a pinned Wigner crystal, the Coulomb repulsion being dominant 

I over the kinetic energy and the disorder. For small r^, the interaction becomes negligible and the 

QQ • ground state is a Fermi glass with Anderson localized one electron states. For intermediate ro. 

CN , a lot of transport measurements following the pioneering works of Kravchenko and co-workers El 

I and made with electron and hole gases give evidence of a metallic phase in two dimensions, 

■ observed □ for instance when 6 < < 9 for a hole gas in GaAs heterostructures. In a recent 
0^ , paper! we have shown that spinless fermions with Coulomb repulsion in disordered 2d clusters 

I exhibits, for a similar range of values, a new ground state characterized by an ordered flow 

' of enhanced persistent currents. A study over the statistical ensemble of the currents supported 

^ ■ by the ground state gives us two well defined values rf and bounding this new phase. 

-a ; 

^ ■ 2 The model 

O ■ 

^ I We consider a simple model of = 4 Coulomb interacting spinless fermions in a random 

■ potential defined on a square lattice with = 36 sites. The Hamiltonian reads: 

<i,j> i 

where cj (cj) creates (destroys) an electron in the site i, t is the strength of the hopping terms 
between nearest neighbours (t = 1 in the following) and is the inter-particle distance for a 
2d torus. The random potential Vi of the site i with occupation number rij = c|cj is taken from 
a box distribution of width W . The interaction strength U yields a Coulomb energy to Fermi 
energy ratio = U / (2ty^7rn7) for a filling factor Tig = N jl? . A Fermi golden rule approximation 
for the elastic scattering time leads to kpl ~ 1927rne{t/W)'^ . Here kp and I denote the Fermi 
wave vector and the elastic mean free path respectively and the above estimate is valid for filling 
factors rie ^ 1. We consider disorder strengths W = 5,10,15 corresponding to kpl = 2.7,0.7 
and 0.3 respectively. The boundary conditions are always taken periodic in the transverse 
y-direction, and such that the system encloses an Aharonov-Bohm fiux in the longitudinal 
x-direction. Imposing (j) = 7r/2 {(j) = n corresponds to anti-periodic boundary conditions), one 
drives a persistent current of total longitudinal and transverse components given by 
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Figure 1: Behavior of a single cluster for W = 15 (left column) and W = 5 (right column) as a function of rs. From 
top to bottom: low energy spectrum (a 1.9rs term has been subtracted); first energy spacing Ao; crystallization 
parameter 7 for the ground state (filled diamonds) and for the first excited state (open diamonds); longitudinal 
Ii (empty circles) and transverse It (empty triangles) current (left scale) and number of occupied sites (filled 

circles, right scale). 



Figure 2; Local currents for the sample of Fig. hi at VK = 5, r-s = (left) and — 17 (right). 



The local current /' flowing at the site i in the longitudinal direction is defined by /' = 
2Im((^'o|c| , , j Cj j l^o))) with an analogous expression for If. The response is paramagnetic 
if > and diamagnetic if Ii < 0. 

3 Numerical results 

Fig. |l| exhibits behaviours characteristic of individual small clusters, comparing the strongly 
disordered regime kpl < 1 {W = 15) with the regime kpl > 1 (W = 5). For kpl < 1, looking at 
the low energy part of the spectrum one can see that, as we gradually turn on the interaction, 
classification of the levels remains invariant up to first avoided crossings, where a Landau theory 
of the Fermi glass is certainly no longer possible. The electronic density pi = (^ol'^il^o) of the 
ground state |\I'o) is maximum in the minima of the site potentials for the Fermi glass. After 
the second avoided crossing, pi is negligible except for four sites forming a lattice of charges 
as close as possible to the Wigner crystal triangular network in the imposed square lattice. 
The degeneracy of the crystal is removed by the disorder, the array being pinned in 4 sites of 
favourable energies. For the same cluster, we have calculated C(r) = J2i PiPi+r and used 
the parameter 7 = max .rC(r) — min.r C(r) to characterize the melting of the crystal (7 = 1 
for a crystal and for a liquid). Calculated for the ground state and the first excited state, 7 
allows us to identify the second crossing with the melting of the crystal. Moreover, one can see 
that the crystal becomes unstable in the intermediate phase, while the ground state is related 
to the first excitation of the crystal. Around the crossings, the longitudinal current /; and 
the participation ratio = ^'^{J2iPi)~^ of the ground state (i.e. the number of sitea^that it 
occupies) are enhanced. The general picture reminds us of strongly disordered chains □ where 
level crossings are associated to a charge reorganization of the ground state and comes with 
an enhancement of the persistent current. For kfl > 1 (W = 5), the previous level crossings 
are now almost suppressed by a stronger level repulsion and charge crystallization occurs more 
continously. There is instead a broad enhancement of /; which, in contrast to the regime kpl < 1, 
is not accompanied by a corresponding increase in the number of occupied sites (,s- One can also 
notice that, increasing r^. It is suppressed before //. This can be understood looking at the local 
currents (see Fig. ^): the interaction causes the transition from a disordered array of currents 
making closed loops to a plastic flow of correlated currents in the longitudinal direction H. 

We summarize in Fig. |3| results from a statistical study of an ensemble of 10^ clusters for 
W = 5, 10, 15. One can see an interaction-induced enhancement of the averaged persistent 




Figure 3: Statistical study of an ensemble of 10^ clusters for W = 5 (circles), 10 (squares) and 15 (triangles) as a 
function of rs- Mean value of the longitudinal current < Ii >. Fraction Cd of diamagnetic samples. Distribution of 
the logarithms of the paramagnetic currents P(Ln|/i,p|) and of the diamagnetic currents P{Ln\Ii^d\) at ra = 1.7 and 
W = 15 (for an ensemble of 10* clusters), fitted by log-normal functions of mean /i — —5.8 and variance — 2.3 
in the paramagnetic CEise, jj, = —6.5 and = 2.9 in the diamagnetic case. Log-averages of the longitudinal 
paramagnetic (empty symbols) and diamagnetic (filled symbols) currents. The straight lines are exponential fits 
giving rY = 9.5 {W — 5), rY = 10.3 iW = 10), rY ~ 12 A {W = 15). Log averages of the transverse currents; 
the exponential fits give = 3.3 {W = 5), rf = 3.7 {W = 10), rf = 4.5 {W = 15). Mean number of occupied 

sites < >■ Mean crystallization parameter < 7 >. 
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Figure 4: Sketch of the phase diagram for 2d spinless fermions in a random potential, with rY (filled circle) and 

(empty circle) values obtained from Fig. ^ 

current < I; > by almost one order of magnitude when ~ 7 and W = b. This effect 
can be partially explained by the suppression of the fraction of diamagnetic clusters The 
paramagnetic Ii^p and diamagnetic Ii^d longitudinal currents, and |/f | have acceptable log-normal 
distributions for all values of Vg when W > f). The stronger is the disorder, the better is the 
log-normal shape of the distributions (see Fig. ^ for W = 15). The averages of the logarithms 
of I;^p, —Ii^d and \It\ give the typical values shown in Fig. ^. They exponentially decay as 
\It\ oc exp(— r^/rf ) and /; oc exp(— r<j/rj^) when is large enough. The values of rf and rY 
extracted from the exponential fits (straight lines of Fig. |3|) are given in Fig. ^ where a sketch 
of the phase diagram is proposed. 

4 Conclusions 

A simple model of Coulomb interacting spinless fermions in a random potential allows us to 
identify a new quantum phase, clearly separated from the Fermi glass and the Wigner crystal. 
This phase is characterized by a plastic flow of currents and a response to an Aharonov-Bohm 
flux with a sign independent of the microscopic disorder. For kpl > 1, displays a strong 
enhancement for intermediate values, which could be the signature of a new metal in the 
thermodynamic limit. For kpl < 1, and // convey similar information while the increase of // 
is accompanied by a decrease of and an increase of 7 when kpl > 1 (see Fig. This suggests 
that, in the latter case, transport in the new phase results from a collective motion of charges 
and not from a delocalization of individual charges. 
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